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Abstract 

The transverse momentum (pj) and centrality dependence of the azimuthal anisotropy (va) are 
measured for charged hadrons in different energies and collision species by the PHENIX exper- 
iment at RHIC. We find that v-i divided by the participant eccentricity of the initial geometry 
proportionally increases with the number of participants to the 1/3 power except at small A^ part 
in Cu+Cu at -^snn = 62.4 GeV. Taking the eccentricity and quark number scaling into account, 
there is a universal scaling for V2 with different energies and collision sizes. The results indicate 
that V2 is determined by more than just the geometrical eccentricity. It also depends on the size 
of collision. We also report that both the behaviors of vi and pj distributions can be understood 
from the thermal nature of produced particles based on hydro-dynamical behavior. 



1. Introduction 

Relativistic heavy ion collisions have been considered as a way to create the quark-gluon 
plasma (QGP), the phase of de-confined quarks and gluons. The Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory was constructed to create and study the QGP. One 
of the most powerful probes for investigating the characteristics of the QGP is measuring the 
azimuthal anisotropy of the charged particles produced by the collisions. The strength of the 
elliptic anisotropy (va) reflects the initial geometrical anisotropy, which creates pressure gradients 
in the collision area. This pressure gradient transfers a spatial anisotropy into a momentum 
anisotropy since particles are pushed harder in the direction of the larger gradients. Thus, the 
measured V2 reflects the dynamical properties of the dense matter produced in the collisions. 

2. Motivation 

One of the most remarkable findings at RHIC is that the strength of V2 can be described 
well by hydro-dynamical models in the low transverse momentum region ( ~ 1 Ge V/c) [ 1 ] . In the 
intermediate transverse momentum region ( 1 ~ 4 Ge V/c), V2 is consistent with quark number (q„) 
and KEj ( = mj - mo) scaling, and the result supports a quark-recombination model [2] . For a 
more comprehensive understanding of v%, we have carried out systematic measurements of V2, by 
measuring x>2 for identified charged hadrons in Au+Au and Cu+Cu collisions at V^nn = 200 and 
62.4 GeV, and studied the dependence on collision energy, species and centrality. We expect that 
the V2 is determined not only by the initial ellipticity, but is also influenced by a finite evolution 
time which can be related to the collision volume. 
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Figure 1: Comparison of integrated vj as a function of A/nan for four collision systems. Left panel shows V2 vs. iVpan, 
middle panel is vj/e vs. Np M and right panel is Vi/(e • ^pan) vs. N pilrt . Closed symbols indicate the results of y^NN = 
200 GeV, open symbols are for y*NN = 62.4 GeV. Circles indicate Au+Au collisions and squares Cu+Cu. The statistical 
and systematic errors are included in the bars. 



3. Results 

3.1. Energy, System Size and Species Dependence 

The left panel in Figure Q] shows V2 vs. A/ part for two collision systems and two collision 
energies. The V2 agrees well at yfs^ = 200 and 62.4 GeV, but there is a clear difference between 
the Cu+Cu and Au+Au results. Since Au+Au and Cu+Cu collisions have different initial geo- 
metrical eccentricities at the same A/ pal1 , applying the eccentricity (e) scaling shows that the data 
follows a universal curve as shown in the middle panel in Figure Q] Here, we use the participant 
eccentricity, calculated with a long and short axis determined by the distribution of participants at 
each collision, using a Monte-Carlo simulation based on a Glauber model and including effects 
from participant fluctuations. The details of the participant eccentricity are described in J3D . One 
can see that v-ils is not a constant and it depends on A/ part . Therefore, vi can be normalized by s 
at the same A/ part , but e is not enough to determine V2 ■ We newly found that v^e is proportional 
to iVi^, and as shown in the right panel in Figure [Tjv2/(e • N^!J t ) is independent of the collision 
system except for small iV part in Cu+Cu at y/s^N = 62.4 GeV. This exception indicates this might 
be a region where the matter has not reached sufficient thermalization, although the errors are 
too large to discuss the difference. A scan of collision energies would be important to further 
study this effect. Figure[TJis for p T = 0.2 - 1.0 GeV/c. The results for p T = 1.0 - 2.0 and 2.0 - 4.0 
GeV/c have the same tendency as well. Similar to the result in Au+Au at ysNN = 200 GeV, V2 
in Au+Au at y/s^N = 62.4 GeV and in Cu+Cu at y^NN = 200 GeV are mostly consistent with 
q„ + KEj scaling for centralities - 50%. In addition to the fact that V2(/?t) does not depend 
on collision energy at RHIC energies, V2 normalized by quark number + KEj, eccentricity, and 
NpJi scaling follows a universal curve as shown in Figure 12 This figure includes the 45 curves 
for 7r/K/p in Au+Au at ys NN = 200 GeV, in Au+Au at \/s NN = 62.4 GeV and in Cu+Cu at 
y/s^ = 200 GeV for the five centrality bins from - 50% in 10% steps. The ^ 2 /NDF of the 
polynomial fitting is 2.1. 

3.2. Radial Flow Effect with Blast-Wave Fit 

To understand the A/ part dependence and KEj scaling behavior for the v'2, we use a Blast-wave 
model to extract dynamical properties of the matter, especially at freeze-out. This model is a two- 
parameter model describing a boosted thermal source based on relativistic hydrodynamics fl. 
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Figure 2: v%/{e ■ N l J^ t ■ q„ ) vs. KEj/q„ for ;r/K/p in Au+Au at -^snn = 200 GeV, in Au+Au at ^Inn = 62.4 GeV and 
in Cu+Cu at y*NN = 200 GeV at five centrality bins for 0-50 % as 10 % step for each system. There are 45 curves. 
Applying polynomial fitting, ^-/NDF is 2. 1 . 



The two parameters, the radial velocity (p\) and the freeze-out temperature (7f ), are extracted 
from the invariant cross section data according to the following equation: 
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where Iq and K\ represent modified Bessel functions with p the transverse boost p(r) = tanh 
=/? s (~), and m T 
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Pj(r) = fisi j), and rrij = ^jp^ + m 2 . Re-plotting the measured pj spectra, weighted by measured 
v 2 (Pt), we obtain the pj spectra in and out-of plane separately for n/K/p. We use two data sets: 
Cu+Cu and Au+Au at v^nn = 200 GeV. Applying the Blast-wave fitting to the pj spectra in and 
out-of plane separately, the fij and 7f in and out-of plane are obtained separately. From these 
values, px2 = (jS T in -pTout)/(j8Tin +pTout)/2 and T fo2 = (T foin - T foout )/(T foin + T foout )/2 are calculated. 
j8t2 (Tfoi) indicates the amplitude of the second harmonic of the /3j (Tf Q ) azimuthal distribution. 
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Figure 3: /?t vs. Np- M for in and out-of plane in Cu+Cu 
and Au+Au at 200 GeV. This result is obtained by fitting 
to PHENIX PRELIMINARY results. 



Figure 4: /?T2/e vs. _/V part in Cu+Cu and Au+Au at 200 
GeV. This result is obtained by fitting to PHENIX PRE- 
LIMINARY results. 
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Figure [3] shows the f3j vs. N pait for in and out-of plane in Cu+Cu and Au+Au collisions. 
The magnitude of j3j is clearly different between in and out-of plane, and it agrees well between 
Au+Au and Cu+Cu, especially for the in-plane data. The amplitude of the second harmonic, /3t2, 
is not the same at the same N pm - t between Cu + Cu and Au + Au. When scaled by eccentricity, 
PtqJs agrees between Au+Au and Cu+Cu and it is fiat at N pirt > 40 as shown in Figure|4]. Since 
V2 is proportional to /?T2 in this model, V2 should be scaled by participant eccentricity. However, 
as shown in Figure[T] this is not what is seen in the measured V2 results. Therefore, this implies 
that v'2 is not determined by only /3x2 but includes other effects such as the freeze-out temperature 
(Tf ) of radial flow. Figure [5] shows the Tf D vs. N pan for both in and out-of plane. The magnitude 
of Tf Q agrees well between Au+Au and Cu+Cu collision at the same N par t. It can be seen that Tf D 
depends on N pmX , and it can influence V2 since the larger Tf D makes pj spectra flatter. Figure [6] 
shows V2 /rig vs. Tf Q obtained by the blast wave calculation, fixing parameters to reasonable 
values «/6t> = 0.5, j8t2 = 0.04, KEj/n q = 0.5 GeV/c) and changing only 7f„. In the freeze out 
temperature region (0.10 < Tf < 0.16 GeV), it can be seen that KEj scaling is approximately 
held by this calculation. In this figure, going from central to peripheral collisions, namely going 
from low Tf to high Tf Q , finite deviations among n/K/p are expected, and V2 for protons becomes 
larger than that for pions. This deviation has the same tendency as the V2 data. 
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Figure 6: V2/n q vs. Tf by the blast-wave calculation with 
Figure 5: r fo vs. N pM for in and out-of plane in Cu+Cu fixed parameters «/? T > = 0.5, /3 T 2 = 0.12, KE T /n q = 0.5 
and Au+Au at 200 GeV. This result is obtained by fitting GeV/c). 
to PHENIX PRELIMINARY results. 
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